Physical ageing studied by a device allowing for rapid thermal equilibration 
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I. INTRODUCTION 
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Ageing of organic glasses to the equilibrium liquid state is studied by measuring the dielectric 
loss utilizing a microregulator where temperature is controlled by means of a Peltier element. Com- 
pared to conventional equipment the new device adds almost two orders of magnitude to the span 
of observable ageing times. Data for five organic glass-forming liquids are presented. The existence 
of an "inner clock" is confirmed by a model-free test showing that the ageing of structure is con- 

■ trolled by the same material time that controls the dielectric properties. At long times relaxation 
' is not stretched, but simple exponential, and there is no "expansion gap" between the limits of the 
, relaxation rates following up and down jumps to the same temperature. 
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■ Ageing is the change of materials properties over time. Ageing phenomena often involve chemical degradation, 
' but there are also many cases of purely physical property changes. A well-known example is the slight shrinking of 
, plastic rulers over years deriving from an extremely slow density increase with time. Understanding physical ageing is 

Q ■ obviously important for applications, but physical ageing also presents fundamental scientific challenges and provides 
c/2 , valuable insight into materials properties. This paper shows that utilizing the Peltier thermoelectric effect ageing 
studies may be extended to considerably shorter times than have so far been possible, adding almost two decades to 
the span of observable ageing times. 

A prime example of physical ageing is that of a viscous liquid's properties relaxing close to its glass transition. In 
equilibrium a liquid's properties do not change with time, of course. If temperature is changed, properties gradually 
^ adjust themselves to new equilibrium values. If temperature is lowered, a glass is produced; recall that by definition a 

Q glass is nothing but a highly viscous liquid that has not yet had time to relax fully to equilibrium [l|, U, H, 0, H, @ . A 
Q , glass gradually approaches the equilibrium liquid state. This state that can only be reached on laboratory time scales, 
however, if the system is kept just below its glass transition temperature ("annealed") - contrary to the popular myth 
I ] windows do not flow observably. 
^ ■ Ageing is an intrinsically nonlinear phenomenon. This is because the ageing rate is structure dependent and itself 
' ages with time 0, H S [O, [13, • Therefore ageing studies provide information beyond those obtained by linear- 
0^ response experiments like dielectric relaxation measurements. A simple ageing experiment consists of a temperature 
step, i.e., a rapid decrease or increase of temperature to a new, constant value. Ideally the new temperature should 
be well-defined homogeneously throughout the sample and constant in time before any relaxation has taken place. If 
that is achieved, it is possible to monitor the complete relaxation to equilibrium of the physical property that is being 
probed. 



II. EXPERIMENTAL 



What are the requirements for such an ideal ageing experiment? First, one should obviously have very good 
' temperature control. Secondly, a physical observable is needed that may be monitored fast and accurately and which 
- preferably - changes significantly for small temperature changes. Third, the set-up should allow for fast temperature 
changes. 

As regards the last point, in current state-of-the-art ageing experiments the characteristic thermal equilibration 
time r is at least 100 s (r is defined by writing the long-time approach to equilibrium oc exp(— t/r)). This reflects 
the fact that heat conduction is a notoriously slow process. Experience shows that in order to monitor a virtually 
complete ageing curve for a temperature down jump, at least four decades of time must be covered. Thus with 
present methods, in the best cases one needs of order 100 s x 10** — 10^ s to perform an ideal temperature down-jump 
experiment. This is more than a week. 

Feynman is famous for his prophesy of the nano-revolution by remarking that there is "plenty of room at the 
bottom." For ageing experiments there is also plenty of "room at the bottom" , but here at the bottom of the time axis: 
Clearly much is to be gained if it were possible to equilibrate sample temperatures much faster. In order to be able to 
make faster temperature-jump experiments we have designed a dielectric cell based on a Peltier thermoelectric element 
by means of which heat flow is controlled via electrical currents (Fig. 1). The characteristic thermal equilibration 
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FIG. 1: Schematic drawing of the dielectric measuring ceh with the microregulator. Liquid samples are deposited in a 50 
fim gap between disks of the dielectric cell. The Peltier element either heats or cools the dielectric cell, depending on the 
direction of the electrical current powering the element. The current is controlled by an analog temperature control system 
that receives temperature feedback information from an NTC thermistor embedded in one disk of the dielectric cell. A stainless 
steel electrode pin keeps the cell pressed against the Peltier element and provides the electrical connection to one of the disks. 



time of the new microregulator is two seconds [14| . This is 50 times faster than the best conventional equipment. If 
required, temperature may be kept constant over weeks, with fluctuations smaller than 100 fiK. 

For monitoring ageing, measuring the dielectric loss at some fixed frequency is probably the best choice. With 
modern equipment this quantity may be monitored fast and accurately. Moreover, for viscous liquids in the right fre- 
quency range the dielectric loss changes considerably for small temperature variations. The use of dielectric relaxation 
measurements for monitoring ageing was pioneered by Leheny et al. p^ . [l6j and by Lunkcnheimer et al. ■ 

Using the microregulator we studied ageing of dibutyl phthalate (DBF), a widely studied glass-forming organic 
liquid with glass transition around 176 K. Ageing was studied by monitoring how the dielectric loss, e"(/ = 0.37Hz), 
develops as a function of time after four 2K temperature steps starting from equilibrium at T = 172, 174, 176 K. In 
order to ensure equilibrium before each temperature jump, the sample was kept at that particular temperature for 
so long that there were eventually no detectable changes of the dielectric properties. Figure 2(a) shows the results, 
jumping from 172 K and 176 K, respectively, to 174 K - as well as the opposite jumps. The inset shows a blow-up of 
the final stages of approach to equilibrium at 174 K. Note that ageing for jumps to a given temperature is considerably 
faster when starting at a higher temperature than when starting at a lower (the two jumps to 174 K, in the middle 
of Fig. 2(a)). This is the so-called fictive-temperature effect described already by Tool in the 1940's |18j . 

Quite generally a relaxation function r[t) is defined by subtracting the long-time limit of the physical property 
monitored. In our case: 



r[t) = |lne"(/ = 0.37Hz,i)-lne"(/ = 0.37Hz,t^cx))|. (1) 
For any relaxation function the Kovacs-McKenna (KM) relaxation rate F is defined 0, El by 

„ d In r 



dt 



(2) 



Figure 2(b) shows the KM relaxation rate plotted as a function of time. At long times there is considerable noise, 
reflecting the fact that the relaxation rate is extremely sensitive to the t ^ oo limit of the dielectric loss being 
determined accurately. Despite this noise it is clear that for the up and down jumps ending at 174 K the KM 
relaxation rates converge to the same number at long times. This shows that there is no so-called expansion gap, as 
Kovacs proposed in 1963 based on monitoring relaxation by measuring refractive index changes Q. The existence 
of an expansion gap was subsequently debated [l^, [13, [2l| , but KoUa and Simon recently concluded that there is no 
expansion gap [221] . The present measurements confirm their simple picture that when equilibrium is approached, 
relaxation rates converge to a unique value. Note that the data of Fig. 2 contradict the widely used stretched- 
exponential relaxation function, exp[— (t/r)^], because according to this function the relaxation rate approaches zero 
at long times, not a constant value. 



III. THE "INNER" CLOCK 



Ageing experiments are usually interpreted utilizing the Tool-Narayanaswami (TN) formalism that models ageing 
in terms of an "inner" clock. This is based on the so-called material time, which is time measured on a clock with a 
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FIG. 2: Ageing of dibutyl phthalate (DBF) studied by monitoring the dielectric loss at 0.37 Hz. (a) Ageing following 2K 
temperature jumps from 172, 174, 176 K. The inset shows the final stages of approach to equilibrium at 174 K from the two 
jumps from 172K and 176 K. (b) The Kovacs-McKenna relaxation rate F defined by eq. ([2]) for the same data. Up and down 
jumps to 174 K at long times give same relaxation rate, which shows that there is no "expansion gap" as proposed by Kovacs 
in 1963 0. (c) Dielectric loss at equilibrium conditions at two higher temperatures where the loss peak is within the measured 
frequency window. This figure illustrates how one utilizes the fact that the loss varies as f~^^^ at high frequencies to measure 
the dielectric loss peak frequency giving the inner clock's clock rate. Details are provided in the supplementary information, 
(d) The KM relaxation rate F (top panel) and its dimensionless version F (bottom panel, defined in the text eq. [S]|n| plotted 
as functions of the normalized relaxation function. The data collapse in the second figure confirms the existence of an inner 
clock (a material time). 



clock rate that changes as sample properties evolve with time [3 HI, 113, HH] • The material time may be thought of 
as analogous to the proper time of relativity theory, the reading on a clock following the observer's world line. The 
TN formalism is standard and widely used in industry for predicting ageing effects. But it is not known whether the 
material clock has real physical existence or is only a convenient mathematical construct. This question cannot be 
answered definitively, but the assumption that an inner clock does exist implies a definite prediction. Thus if a material 
time exists controlling all relaxations, the relaxation of the clock rate activation energy must be determined by the 
same material time that controls the dielectric ageing process. We proceed to formulate a quantitative, model-free 
test of this. 

The dielectric loss of DBP has a high-frequency power-law decay described by e"(/) oc j^^l"^ [2^. [27|. Suppose the 
clock rate ^(t) controls both structural relaxation and dielectric relaxation. We normalize 7(t) such that it equals the 
loss peak frequency when the liquid is in equilibrium. Details are given in the Supplementary Information, but the 
main idea is illustrated in Fig. 2(c): Even out of equilibrium the dielectric relaxation rate ^(t) may be determined 
by measuring the high-frequency loss utilizing 



."(/,t) « (^)'^'. (3) 

If structural relaxation is also controlled by 7(t), this expression allows us to determine the structural relaxation clock 
rate via 
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lnj{t) = 2lne"{f,t)+A. (4) 

The constant A is found by cahbrating using cquihbriuni data from higher temperatures where the loss peak is within 
the observable frequency range (Fig. 2(c)). 

The material time i is the actual time measured in units of the (time-dependent) inverse clock rate [l^, [l^, [13, [lH, 
i3, [11,111,111, di=dt/j{t)~\ i.e. 



di = j{t)dt. (5) 

According to the TN formalism, for all temperature jumps the normalized relaxation function R{t) = r{t)/r{0) is 
some function of the material time that has passed since the jump, a function that is independent of the sign and size 
of the temperature jump: R = R{t). Generally the function R{t) depends on the quantity that is monitored. 

The existence of an inner clock, i.e., the assumption that the dielectric clock rate equals the structural relaxation 
clock rate, can be checked without evaluating t explicitly or fitting data with analytical functions. This is done by 
proceeding as follows. First, define the dimensionless KM relaxation rate by 



7(*) 

This quantity is given directly by experiment since both T{t) and j{t) are determined from lne"(/, t). Next, we note 
that r = —d\nr/dt = —dlnR/dt. This implies that if R{t) is a unique function of t, i.e., independent of sign and size 
of the temperature jump, then F = T{t) is also unique. In turn this implies that F is a unique fmiction of R: 



f = ^R) . (7) 

Figure 2(d) shows that this prediction is fulfilled within the experimental uncertainty: The upper part of that figure 
plots the KM relaxation rates F as functions of R for the four temperature jump experiments of Fig. 2(a), the lower 
part plots F as functions of R. The long time limit corresponds to i? ^ 0. The fact there is data collapse confirms 
the existence of a material time that controls both dielectric ageing and the ageing of the relaxation rate itself (the 
Supplementary Information provides mathematical details). 

We repeated the experiment for four other liquids, diethyl phthalate, 5 poly-phenyl ether, 2,3 epoxy phenyl propy- 
lether and triphenyl phosphite. These liquids all have high-frequency dielectric losses proportional to /"^^^ |26l. [27j. 
The results for the KM rates F(i?) are plotted in Fig. 3(a) and the corresponding dimensionless KM rates in Fig. 
3(b). The spread in KM rates at i? = reflects the already mentioned fact that relaxation rates are almost impossible 
to determine accurately at long times where the noise is of the same magnitude as the distance to equilibrium. In all 
cases we find good data collapse, confirming the existence of a material time for these liquids. 



IV. DETAILS OF THE TEST 



We here detail the assumptions that must be made in order to arrive at the proposed model-free test of the 
material clock hypothesis. We first define a generally time-dependent dielectric relaxation rate 7d(i) by generalizing 
time-temperature superposition, then define the structural relaxation rate 7s (Oi finally arrive at the test. 



A. The dielectric relaxation rate defined for out-of-equilibrium situations 

The simple power-law high-frequency dielectric loss of DBF and the other liquids studied here, e"(/) oc /^^^^ 
[IE [13 1 is utilized to monitor the dielectric relaxation rate jdit) as the structure ages following a temperature change 
by proceeding as follows. By definition 7^ = fp in the equilibrium liquid phase where fp is the dielectric loss-peak 
frequency. If temperature is lowered in a step experiment, the dielectric loss curve continuously moves to lower 
frequencies as the system relaxes to equilibrium. The idea is now (compare Fig. 2(c)) that how much the dielectric 
relaxation rate 7^ has changed, i.e., how much the loss peak is displaced, may be determined from how much the loss 
at some fixed frequency in the high-frequency power-law region changes. Thus if the loss is continuously probed at a 
fixed frequency, the dielectric relaxation rate is determined from 
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FIG. 3: Results for four other liquids, (a) Kovacs-McKenna relaxation rates as functions of the normalized relaxation function 
for the following four liquids: Di-ethyl phthalate (DEP); 5-poly phenyl ether (5PPE); 2,3 epoxyphenyl-propylether (2,3 epoxy); 
tri- phenyl- phosphite (TPP). (b) Dimensionless Kovacs-McKenna relaxation rates as functions of the normalized relaxation 
function for the same liquids. In all cases there is good data collapse within experimental errors, confirming the existence of 
an inner clock (a material time) for these liquids. Note that for TPP even the 4K down-jump follows the same master-curve 
as the 2K up- and down jumps. For all four liquids the dielectric loss was monitored at frequency f — 1 Hz. 
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e"if,t) = aif/Mt))-'^^- (8) 

Although this idea seems fairly straightforward, the concept of a dielectric relaxation rate has no obvious definition in 
a non-stationary situation as that of a relaxing structure. In order to specify the precise assumption needed to define 
7d(t) and justify eq. ([8]), we reason as follows: According to standard linear-response theory, for a system in thermal 
equilibrium the measured output is calculated from a convolution integral involving the history of the input previous 
to the measuring time. A convenient way to summarize time-temperature superposition for the equilibrium liquid 
is to formulate the convolution integral in terms of a dielectric "material" time t: If 7^ is the equilibrium liquid's 
dielectric relaxation rate (7^ — fp), the dielectric material time is defined from the actual time thyt = ^^t. In terms 
of t, if the input variable is the electric field E and the output is the displacement vector D, the convolution integral 
is of the form 



D(i) = / -E{i-i')e{i')di' . (9) 
"'0 

Equation ^ describes time-temperature superposition because it implies that, except for an overall time/frequency 
scaling, the same frequency-dependent dielectric constant is observed at different temperatures (we ignore the tem- 
perature dependence of the overall loss, an approximation which introduces an error of just 1% over the range of 
temperatures studied). 

In eq. © the dielectric material time is defined from the actual time by scaling with 7^, a quantity that is of course 
strongly temperature dependent. In the out-of-equilibrium situations following a temperature jump, the simplest 
assumption is that eq. ([9|) also applies during the ageing process, but with the dielectric material time generalized by 
now assuming a time-dependent dielectric relaxation clock rate ^dit): 



dt^-id{t)dt. (10) 

Ks t 00 7rf(t) approaches the equilibrium liquid's loss peak frequency fp at the new temperature, of course. 

The equilibrium liquid's power-law dielectric loss e" oc /^^^^ applies whenever / ^ fp. Equation ([9|) thus implies 
that tit') oc whenever t' <C 1. It now follows that the proposed generalization of eq. ([9]) now implies that the 

dielectric relaxation rate ^d{t) may be determined from of eq. ([5]). 



B. The structural relaxation rate 



We define the dielectric relaxation rate activation (free) energy E, which generally depends on time, by writing 
(where /3 = l/ksT and the prefactor 70 is taken to be 10"^'' s) 



7<i(t) = 7oexp(-/3S(i)). (11) 

The activation energy depends on structure and it relaxes during strucural relaxation. From the Tool-Narayanaswami 
formalism one expects that after temperature changes the activation energy relaxes following a convolution inte- 
gral over the temperature history, with a properly defined material time. For convenience we include the inverse 
temperature in the defining equation by writing for small temperature variations AT{t) 



A{f3E){i) ^ / AT{t-t')(j){t')dt' . (12) 
Jo 

Here the structural relaxation material time is defined by 



di^-1s{t)dt, (13) 

where 7s {t) is the rate of structural relaxation (one may have different structural relaxation rates and thus different 
material times for differing physical quantities relaxing (volume, enthalpy, ...)). 
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C. Assuming the existence of a material cloclc 



The assumption that the material clock has real physical existence implies that the dielectric relaxation rate is 
identical to the structural relaxation rate: 



ls{t) = jd{t) ■ (14) 

This is not a trivial assumption. Thus Eqs. © and (fT^ may both well apply, but with different definitions of the 
reduced time. If Eq. (fT4|) is obeyed, however, Eq. ([8]) implies that after a temperature jump the logarithm of the 
measured loss, 



ln6"(/,t) = ^+C, (15) 

relaxes with the rate that is determined from itself via Eq. ([8]). It is this prediction that is tested in the main paper's 
Eq. (O: The dimcnsionless KM relaxation rate F of Eq. © is defined by dividing the KM relaxation rate F by 
the structural relaxation rate 7s; this quantity however is determined by means of Eq. ^ that actually gives the 
dielectric relaxation rate 7^. 



V. SUMMARY 



By measuring several relaxing quantities the existence of a material clock can easily be checked by simply investi- 
gating whether or not the quantities relax following the same function of time for various temperature jumps. This is 
a direct, model- free test, i.e., one that neither involves free parameters nor the fitting of data to some mathematical 
expression. For measurements of just one relaxing quantity as in this paper we believe that the above is the first 
model-free demonstration of the existence of an inner clock. 

By modern micro engineering it should be possible to extend ageing experiment to even shorter times, thus making 
it realistic to perform a series of ideal temperature-jump experiments over just hours. Thus it is not unlikely that - 
eventually - ageing studies could become routine on par with, e.g., present-day dielectric measurements. 
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